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ABSTRACT 

Irrigation is an important human activity that may impact local and regional climate, but current climate 
model simulations and data assimilation systems generally do not explicitly include it. The European Centre 
for Medium-Range Weather Forecasts (ECMWF) Interim Re- Analysis (ERA-Interim) shows more irriga- 
tion signal in surface evapotranspiration (ET) than the Modern-Era Retrospective Analysis for Research and 
Applications (MERRA) because ERA-Interim adjusts soil moisture according to the observed surface 
temperature and humidity while MERRA has no explicit consideration of irrigation at the surface. But, when 
compared with the results from a hydrological model with detailed considerations of agriculture, the ET from 
both reanalyses show large deficiencies in capturing the impact of irrigation. Here, a back- trajectory method is 
used to estimate the contribution of irrigation to precipitation over local and surrounding regions, using 
MERRA with observation-based corrections and added irrigation-caused ET increase from the hydrological 
model. Results show substantial contributions of irrigation to precipitation over heavily irrigated regions in 
Asia, but the precipitation increase is much less than the ET increase over most areas, indicating that irri- 
gation could lead to water deficits over these regions. For the same increase in ET, precipitation increases are 
larger over wetter areas where convection is more easily triggered, but the percentage increase in pre- 
cipitation is similar for different areas. There are substantial regional differences in the patterns of irrigation 
impact, but, for all the studied regions, the highest percentage contribution to precipitation is over local land. 
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1. Introduction 

Irrigation is an important human activity that has the 
potential to impact local and regional climate through 
the hydrological cycle and surface energy balance (e.g., 
Chase et al. 1999; Pielke et al. 2011). About two-thirds 
of the global freshwater withdrawals from surface and 
underground are used for agriculture (Shiklomanov 2000), 
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of which irrigation is the main component. Although 
there are still uncertainties in the amount of irrigation 
over various regions (Wisser et al. 2008), numerous 
modeling and observational studies have consistently 
shown that irrigation can change the surface energy 
partition by reducing sensible heat flux and increasing 
latent heat flux, thereby reducing surface temperature 
and the diurnal temperature range (e.g„ Adegoke et al. 
2003; Haddeland et al. 2006; Mahmood et al. 2006; 
Kueppers et al. 2007; Sacks et al. 2009; Lobell and 
Bonfils 2008; Ozdogan et al. 2010; Sorooshian et al. 2011). 
In addition to the impact on surface climate, irrigation 
also increases the amount of water vapor in the atmo- 
sphere through evapotranspiration (ET) (Boucher et al. 
2004; Mahmood et al. 2008). The additional moisture 
can enhance the convective available potential energy 
(CAPE) within the boundary layer and make 
it thermodynamically more conducive to increase in 
rainfall (Pielke and Zeng 1989; Betts et al. 1996). These 
have been demonstrated by some modeling (e.g., Segal 
et al. 1998) and observational (e.g., Moore and Rojstaczer 
2002; DeAngelis et al. 2010) studies. Negative impacts 
of irrigation on precipitation, which have also been 
found in some modeling studies, are mostly caused by 
a weakened monsoon circulation (e.g., Saeed et al. 
2009; Douglas et al. 2009; Puma and Cook 2010; 
Guimberteau et al. 2012), although negative soil 
moisture-precipitation feedback is likely in some 
conditions (Findell and Eltahir 2003). 

Therefore, irrigation can contribute to precipitation 
directly as precipitable water vapor and at the same time 
change the atmospheric instability, or it can impact 
precipitation indirectly by modifying the atmospheric 
circulation and moisture transport. These effects usually 
happen concurrently during large-scale irrigations (e.g., 
Saeed et al. 2009; Douglas et al. 2009; Lee et al. 2011). 
Observational studies cannot easily separate them, 
while model simulations, even if capable, are not com- 
pletely reliable. As indicated by some land model in- 
tercomparison experiments (e.g., Henderson-Sellers 
et al. 1996; Dirmeyer et al. 2006), the response of ET 
to the availability of surface water is highly model 
dependent. Coupling land models with atmospheric 
models can amplify the uncertainties (Wei et al. 2010). 

in this study, we use reanalysis data and estimated ET 
over irrigated areas to examine the direct impact of ir- 
rigation on precipitation, that is, through its contribution 
to precipitable water vapor and the associated impact on 
convective instability. A back-trajectory method is used 
to investigate the water vapor transport. Irrigation may 
substantially modify the large-scale circulations and 
patterns of water vapor transport, but the reanalysis data 
should implicitly include these changes as a result of the 


assimilation of meteorological observations. Therefore, 
we estimate the impact of irrigation on precipitation 
under these already-modified climate conditions, which 
is different from the total impact of irrigation as assessed 
by some model simulations (e.g., Saeed et al. 2009; 
Douglas et al. 2009; Lee et al. 2011; Guimberteau et al. 
2012). This separates the direct impact of irrigation on 
precipitation from other indirect impacts. 

2. Data and method 

a. Irrigation area map 

The Food and Agriculture Organization of the United 
Nations (FAO)/University of Frankfurt global map 
of irrigated areas (Siebert et al. 2005, 2007) is used 
to define irrigation area (Fig. 1). It reports the per- 
centage of area equipped for irrigation around the year 
2000. The original data has a resolution of 5' (V 12 0 ), and 
Fig. 1 shows the data at V 2 0 resolution. Four intensively 
irrigated regions — northern India/Pakistan, the North 
China Plain, the California Central Valley, and the Nile 
Valley — are selected for specific analysis. 

b. Reanalysis data 

In this study, the Modern-Era Retrospective Analysis 
for Research and Applications (MERRA; Rienecker 
et al. 2011) is primarily used. MERRA was generated 
with the Goddard Earth Observing System (GEOS-5) 
atmospheric model and data assimilation system, with 
a particular emphasis on the hydrologic cycle. For the 
land variables, we use an early version of a product called 
“MERRA-Land” (Reichle et al. 2011) instead of the 
original MERRA data. MERRA-Land is more accurate 
than the original MERRA due to a correction to the 
model-generated MERRA precipitation with the 
observation-based Global Precipitation Climatology 
Project (GPCP) pentad product (Huffman et al. 2009; 
Xie et al. 2003) and revision of parameter values in the 
original canopy precipitation interception model, which 
caused large biases. MERRA and MERRA-Land do 
not consider irrigation in the land model and also do not 
assimilate surface observations of temperature and 
humidity. MERRA covers the period from 1979 to the 
present and MERRA-Land covers the period from 1979 
to 2008. The horizontal grid of both is 2 /3° longitude by 
Vi latitude. 

The European Centre for Medium-Range Weather 
Forecasts Interim Re-Analysis (ERA-Interim; Dee 
et al. 2011) is used for comparison with MERRA. 
Although ERA-Interim also does not explicitly include 
irrigation, it uses surface observations of temperature 
and humidity to update soil moisture and soil tempera- 
ture estimates for each of the four layers of the land 
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Fig. 1. The percentage of area equipped for irrigation around year 2000 from the FAO/ 
University of Frankfurt (Siebert et al. 2005, 2007). The four regions in the red boxes are for 
detailed analysis. 


surface model (Douville et al. 2000; Mahfouf et al. 
2000). This indirectly considers the effect of irrigation 
on soil moisture and ET (Tuinenburg et al. 2012). The 
ERA-Interim data cover the period from 1979 to 
2010 and have a horizontal grid of 0.75° X 0.75°. Both 
MERRA and ERA-Interim are examples of a new 
generation of reanalysis that incorporates many im- 
portant model improvements and changes in analysis 
methodology. 

Figure 2 compares the ET from MERRA-Land and 
ERA-Interim over four selected irrigation regions. Each 
region (land points only and excluding Tibet) is separated 
into two halves with low and high irrigation based on the 
percentage of irrigated area in each grid box (shown in 
Fig. 1). ET over the two areas is shown separately to 
identify the impact of irrigation. For both MERRA and 
ERA-Interim, the more irrigated areas generally show 
higher ET than the less irrigated areas, even though there 
is no explicit irrigation in MERRA. (Note that there may 
be natural climate differences between the two separated 
areas in each region that contribute to the differences in 
ET, and the assimilation of other observational data may 
also have some impact on ET.) However, the ET dif- 
ferences between the two separate areas are much 
larger in ERA-Interim than in MERRA, especially 
over northern India and the Nile Valley, showing the 
evident irrigation signals. These differences could be 
caused by the assimilation of observed temperature 
and humidity, such that observed cooler temperature 
or higher humidity over irrigated areas is compensated 
by added soil moisture in ERA-Interim (Tuinenburg 
et al. 2012). 

c. A simple estimate of irrigation-caused ET increase 

As there is no explicit irrigation over land in MERRA 
or MERRA-Land, we assume that the irrigation signals 


in atmospheric forcings have little impact on ET 
(shown below) and test a simple method to estimate 
the impact of irrigation on ET. For each land grid 
point, the total land area includes an irrigated portion 
/ and a nonirrigated portion 1 — /. The total evapo- 
transpiration (Etotai) then comes from both of the land 
areas: 


^total ^irrigated^ 


-irrigated^ /)’ CO 


where / can be obtained from the irrigated area data- 
set, the ET rate over nonirrigated land (£non-imgated) is 
simply the MERRA-Land output, and the ET rate over 
irrigated land is estimated as 


^irrigated 


= (SH + LH)EF, 


( 2 ) 


where SH + LH is the sum of latent plus sensible heat 
fluxes and is from MERRA-Land and EF is evaporative 
fraction (the ratio of LH to SH + LH). For irrigated land 
EF is specified as 0.8, an estimated value from model 
simulations (Dirmeyer et al. 2000). Therefore, the ET 
increase due to irrigation (dE) is 


^irrigation ^'irrigated ^non-irrigated 


)/,0]. (3) 


Evidently, 


F = F + F 

total non-irrigated irrigation * 


( 4 ) 


Figure 3 shows the estimated monthly ET climatology 
over the four selected irrigation regions for irrigated 
(£)otai) an d nonirrigated (E n0 n-irrigated) conditions and 
their difference (Eirrigation)- The highest estimated 
irrigation generally happens in the dry and warm 
season, when the available energy for ET is high but 
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Fig. 2. Monthly ET climatology from (left) MERRA-Land (1979-2008) and (right) ERA-Interim (1979-2010) over two separate halves 
with high and low irrigation for each region defined in Fig. 1. Ocean grid points and Tibet (altitudes above than 3000 m) were removed 
before separating the areas. 


the soil does not contain much water, crudely de- 
scribing the role of irrigation to provide water for 
agriculture when demand exceeds the rainfed supply. 
When the rainy or cold season begins, the irrigation 
drops. 


d. ET from a hydrological model with consideration 
of human activities 

The method above gives a simple estimate of dE 
based on energy constraints. However, such a simple 
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California Nile Valley 




Fig. 3. Estimated monthly ET climatology with and without irrigation and their difference (ET increase due to irrigation) during 1979- 
2008 from a simple method and based on MERRA-Land data (see section 2c). The results are only for the more irrigated halves of the four 
regions (as defined for Fig. 2). 


estimate neglects the surface water requirements 
of different crop types and crop production and also 
neglects the temporal variations in crop irrigation de- 
mand as a result of planting and harvesting. In this 
study, the ET data from a global reconstruction of 
hydrography by Wisser et al. (2010) are used to eval- 
uate the simple estimate. The reconstruction uses 
WBMplus (Vorosmarty et al. 1998), a water balance 
and transport model that simulates daily irrigation 
water demand by combining data on irrigation area, 
growing seasons, cropping patterns, soil properties, and 
climate. The simulated mean annual irrigation water 
demand, aggregated by country, correlates well with 
national statistics reported by the FAO (Wisser et al. 
2008). Because of these sophisticated considerations 
and some proven success, the ET data from the 
hydrological model are used for further analysis. 

Monthly ET data for 1979-2002 from two simulations 
with the model, a preindustrial run (without irrigation) 
and a contemporary run (with irrigation), are obtained. 
Figure 4 shows the reconstructed climatological ET 
from these two simulations and their difference dE. 
Compared to the simple estimate shown in Fig. 3, the ET 
without irrigation is similar, but dE is very different, 
much larger in amplitude, and shifts in seasonality, 


except over California where the growing season is 
long. Wisser et al. (2008) found that the irrigation 
requirement is very sensitive to the factors related to 
paddy rice. This may be the reason for the high dE 
over northern India during autumn when the autumn/ 
winter rice is planted, and this could also be the reason 
for very high dE over the Nile Valley. The higher dE 
over northern China during spring and early summer 
could be caused by the low irrigation efficiency; that 
is, a large amount of water is lost during conveyance 
and much of it is transformed to ET (Doll and Siebert 
2002 ). 

By comparing Fig. 4 with Fig. 2, we find that the 
impact of irrigation on ET is very little captured by the 
reanalyses, especially MERRA, although the atmo- 
spheric forcings may contain some irrigation signals 
from data assimilation. This should be mainly caused by 
the lack of sophisticated agricultural modeling and 
land assimilation (e.g., for rice paddy) in the reanalyses, 
although the estimated irrigation withdrawal from this 
hydrological model is a little higher than that from other 
studies (Guimberteau et al. 2012). We can therefore 
assume that the ET from MERRA-Land contains no 
irrigation, and the irrigation-caused ET increase can be 
added to it to represent the effect of irrigation. 
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Fig. 4. As in Fig. 3, but from the preindustrial (ET without irrigation) and contemporary (ET with irrigation) simulations of 

Wisser et al. (2010) for 1979-2002. 


e. Quasi-isentropic back-trajectory method 
and data used 

The quasi-isentropic back-trajectory (QIBT) method 
(Dirmeyer and Brubaker 1999, 2007; Brubaker et al. 
2001) is used to track the water vapor for each pre- 
cipitation event backward in time along the isentropic 
surfaces, assuming precipitated water is drawn from 
the atmospheric column in a distribution that follows 
the vertical profile of specific humidity. Traces start 
from the grid box that has precipitation and from the 
time step when precipitation occurs, backward in 
space and time until all of its original precipitation is 
attributed to ET, but no longer than 15 days (the av- 
erage residence time of moisture in the atmosphere is 
about 10 days). The time step for calculation is 45 min, 
and the calculations and output are aggregated into 
pentads. See the reference papers above for a detailed 
description of the QIBT method. An advantage of the 
QIBT method compared to using water vapor tracers 
in models (e.g., Druyan and Koster 1989; Numaguti 
1999; Bosilovich and Schubert 2002) is that it can be 
applied a posteriori to 3D analysis data, and the study 
region can be defined after the model simulation has 
been finished. Several similar methods have been 
developed to examine the water vapor sources for 


precipitation (e.g., Stohl and James 2004; Sodemann 
et al. 2008; van der Ent et al. 2010). 

The required data for the QIBT calculation include 
precipitation and ET at the surface and temperature, 
humidity, and wind at different tropospheric levels, 
preferably at a subdiurnal time step (6-hourly data are 
used in this study). All data are from the original 
MERRA except ET and precipitation. Since we are 
concerned with the impact of irrigation, monthly dE 
from the hydrological model (section 2d) is added to 
the daily ET data from MERRA-Land by scaling, and 
this new ET data are used to correct the 6-hourly 
MERRA ET. Daily gauge-based Climate Prediction 
Center (CPC) unified precipitation (Xie et al. 2007; 
Chen et al. 2008) is used to correct the MERRA 
precipitation. It has a spatial resolution of 0.5° X 0.5° 
and covers global land. The two corrections are made 
at the pentad time scale by multiplying the original 
MERRA variable at each time step with a ratio of 
the pentad-average value of the target (MERRA- 
Land ET plus irrigation or CPC unified precipitation) 
and the pentad-average value from original MERRA. 
Note that small negative ET changes also happen over 
some areas. They are considered in the QIBT run, but 
ET from these areas is not counted as a contribution 
from irrigation. 
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a) Annual ET increase (mm) 



b) Percentage of ET increase (%) 



Fig. 5. (a) Ten-year (1986-95) mean annual total ET increases due to irrigation (the dif- 
ference between the contemporary and preindustrial simulations of Wisser et al. 2010). (b) The 
percentage of ET increases relative to the mean preindustrial condition. Almost identical re- 
sults are found for the 24-yr (1979-2002) average. 


The QIBT run is performed for the 10-yr period 1986- 
95. Calculations are only made for the northern low- and 
midlatitude (5°-55°N) land where most irrigated areas 
are located. The output is a spatial distribution of the 
evaporative sources for precipitation over each grid 
point (cf. Dirmeyer and Brubaker 1999, 2007; Dirmeyer 
and Kinter 2010; Wei et al. 2012). We can calculate the 
impact of dE over a certain region on precipitation over 
local and surrounding regions by aggregating the evap- 
orative contribution from this region. In estimating the 
impact of dE , we assume dE accounts for a fixed portion 
of the total ET over each grid point during a month. 

3. Results 

Figure 5 shows the mean annual dE, which is the dif- 
ference between the contemporary run and preindus- 
trial run of the hydrological model, as introduced in 


section 2d. For the convenience of comparison, data for 
the same 10-yr period as the QIBT run is shown. The 
spatial pattern of dE is very similar to the percentage of 
land dedicated to irrigation shown in Fig. 1. Figure 5b 
demonstrates that after irrigation is considered ET is 
almost doubled over northern China and northern India 
and more than tripled over northwestern India and 
Pakistan, the Nile Valley, and some small areas over 
the Middle East, California, etc. This illustrates the 
importance of human activities on surface ET over some 
regions, although the global average dE over land 
(shown in Fig. 5) is only about 2.5%. 

Figure 6 shows the amount of precipitation coming 
from irrigation (dP) and its fraction in total precipitation 
(dP/P) estimated from the 10-yr QIBT analysis. The dP 
is high over India and northern Pakistan, central and 
northern China, and northwestern Ethiopia, basically 
consistent with the pattern of ET increase but also 




282 


JOURNAL OF HYDROMETEOROLOGY 


Volume 14 


a) 


Annual precip. from irrigation (mm) 



b) 
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Percentage of total precip. (%) 



1 2 3 4 5 7 10 15 20 25 

Fig. 6. Results of the 10-yr QIBT analysis that shows (a) the mean annual total precipitation from irrigation and 
(b) the percentage of the total precipitation that is from irrigation. 


representing the nonlocal moisture transfers. The fraction 
of total precipitation that is dP shows a different pattern, 
with high values over northwestern India and Pakistan, 
western China, and northern China. On average, only 
1.8% of the land precipitation (5°-55°N) is from irri- 
gation. The average dE over northern land (5°-55°N) 
is 22.71 mm yr -1 (shown in Fig. 5a), but dP over the 
same area is 10 mm yr -1 . As the long-term moisture 
change in the atmosphere is negligible, the remainder of 
dE should mostly contribute to the precipitation over 
ocean. 

In addition to dE , the patterns in Fig. 6 are related to 
other factors such as the wind field and climate over 
local and surrounding regions. The strength and pattern 
of the wind field determines whether the impact is local 
or more remote, and the regional climate conditions 
control the impact of ET increase on precipitation 
(Findell et al. 2011; Wei and Dirmeyer 2012). To ex- 
amine these impacts, we investigate the relationship 
between dP and dE. Figure 7a shows the ratio of dPIdE 
during March-November; the boreal winter months are 
excluded owing to the weak impact of ET on precip- 
itation. The data are interpolated to 5° X 5° grids to 
highlight regional budgets. The irrigation regions with 
dPIdE less than 1 (precipitation increase is less than the 
water loss) may experience water deficits if there are 
not enough renewable water resources (e.g., rivers), 
while regions with dPIdE > 1 are less likely to have 
water deficits. It can be seen that dPIdE is less than 1 
over most irrigated areas, consistent with the results 
above. The heavily irrigated regions usually have a small 


value of dPldE\ dPIdE is greater than 1 over a few 
regions, including southeast India, southwest China, 
northeast China, and the Korean peninsula. They are 
usually downwind of the regions with high dE but have 
small dE themselves, indicating that they are not heavily 
irrigated. According to Wada et al. (2010), these regions 
are a few that are not under large groundwater de- 
pletion, supporting our results. 

Separating the irrigated grids in Fig. 7a into dry and 
wet halves, Fig. 7b shows that dP increases with dE , 
and for the same dE wet areas have significantly larger 
dP than dry areas. Figure 7c further shows that the 
percentage precipitation increases ( dP/P ) are actually 
very similar over dry and wet areas, and they both in- 
crease with dE. This means that the relative contribution 
of dE to precipitation is similar over all regions and 
higher dE leads to higher dPIP. Evapotranspiration 
increases over dry regions struggle to produce enough 
CAPE to trigger local precipitation, and the moisture is 
instead advected away. In humid regions the ET 
increase can more easily enhance convective instability 
and contribute to increased precipitation (e.g., Dirmeyer 
et al. 2010). Our results are supported by those from 
some numerical experiments. For example, Segal et al. 
(1998) found that the precipitation increase in response 
to irrigation typically occurs over areas with existing 
precipitation, and Ter Maat et al. (2006) showed that the 
large-scale irrigation over southwest Saudi Arabia 
can lead to very limited water recycling and rainfall 
generation. Note that we emphasize the impact of ET 
increase on precipitation in the perspective of water 
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Fig. 7. (a) Ratio of March-November total precipitation increase ( dP ) and ET increase (dE) due to irrigation over 
5° X 5° grid boxes. Only grid boxes with dE > 10 mm are shown, (b) Scatterplot of dP and dE for each grid box in (a), 
(c) Scatterplot of dPIP and dE for each grid box in (a). The grid boxes are separated into two equal-sized groups (red 
and blue colors) according to their total precipitation, and their respective linear regressions are shown (they are not 
straight in the figure because of the log scale). The dashed lines are their respective 95% confidence intervals. The 
regressions are significant at p < 0.01 level even when using a conservative estimate of the degree of freedom of one- 
third of the sample size due to spatial correlations. 


cycle, that is, by contribution of water vapor, which is 
somewhat different from the impact of ET on the variation 
and prediction of precipitation (Wei and Dirmeyer 2012). 

Next, the contributions of irrigation to precipitation 
are calculated and discussed separately for each of the 
four heavily irrigated regions. Figure 8 shows the con- 
tribution of irrigation over northern India to regional 
precipitation. This region is the most heavily irrigated 
region in the world as measured by both the percentage 
of land equipped for irrigation (Fig. 1) and dE (Fig. 5). 
Satellite-based observations have shown strong 
depletion of groundwater over this region in recent 
years (Rodell et al. 2009). Figure 8 shows that, although 
irrigation has a strong impact on local precipitation, it 
also affects precipitation outside of the irrigation region, 
especially downwind (Figs. 8a, b). The biggest impact on 
dP is during August-September (Figs. 8c,d), although 
dE is largest during September-October (Fig. 4). This is 
because August is much more humid than October and, 


according to the mechanism discussed above, it is much 
easier for enhanced precipitation to be triggered. The 
dPIP is highest during October (Figs. 8e,f) with a local 
(inner box) average of 15.2% and high values above 
60% over northwest India and eastern Pakistan (not 
shown), but the actual precipitation is low during this 
time (less than 5 mm month -1 ). 

The north China plain is one of the main agricultural 
regions in China and has a very dense population. Irri- 
gation over this region also affects local and remote pre- 
cipitation but is smaller in amplitude than over northern 
India (Fig. 9). This is related to both the irrigated area 
and irrigation amount, which are both smaller. A large 
amount of the irrigation over this region is transported 
to the northeast by the monsoon flow during June-July 
when irrigation and its impacts are largest. The dPIP is 
highest during June with a local average of 6.3% and 
high values above 12% over southeast Beijing, north- 
west Tianjin, and part of Hebei province. 



284 


JOURNAL OF HYDROMETEOROLOGY 


Volume 14 


a) Annual precip. from irrigation (mm) b) Percentage of total precip. (%) 



1 3 5 7 9 11 13 14 
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2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 

Fig. 8. (a) The results from 10-yr QIBT analysis showing the mean annual total precipitation (mm) from northern 
India irrigation in the inner red box (shaded). The outer red box is the area for zonal or meridional averages in (c)-(f). 
The arrows show the annual average vertically integrated moisture flux (kg m -1 s -1 ). (b) Percentage of the total 
precipitation that is from irrigation in the inner red box. (c) Latitude-time plot for the monthly climatology of zonal 
average contribution of irrigation to precipitation (mm month -1 ). The horizontal dashed lines are the north and 
south boundaries of the source region [inner red box in (a)], (d) As in (c), but the longitude-time plot for the 
meridional averages. The vertical dashed lines are the west and east boundaries of the source region, (e) As in (c), but 
shown as the percentage of total precipitation, (f) As in (e), but the longitude-time plot. 


The California Central Valley is a major agricultural 
region in the United States. Irrigation water, mainly 
from snowmelt and groundwater, is very important for 
agriculture in this region. Figure 10 shows the contri- 
bution of irrigation over this region to precipitation. 
It can be seen that the contributions are much weaker 


than for northern India or northern China, which is 
related to the smaller irrigated area and the dry climate 
over surrounding regions. The location that is most 
affected in amount moves from local areas in May to 
remote locations (Montana, Wyoming, South Dakota, 
and North Dakota) by July. The local dPIP is highest 


February 2013 


WEI ET AL. 


285 


a) Annual precip. from irrigation (mm) 



c) Zonal average (mm/mon) 
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b) Percentage of total precip. (%) 
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d) Meridional average (mm/mon) 
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1 1.5 2 2.5 3 3.5 4 4.5 5 


Fig. 9. As in Fig. 8, but for irrigation over northern China. 


during July with an area average of 1.9% and peak 
values over 10%. 

The Nile Valley is a unique agricultural region that 
lies within a severe desert. It mainly depends on water 
from the Nile River for irrigation. Because of the very 
dry climate over this region, the evaporated irrigation 
water is transported far afield, south to Ethiopia and 
southern Sudan and northeast to western Asia, where 
the climates are more humid (Fig. 11). The most affected 
region is Ethiopia during July-August, but it only 
accounts for a very small percentage of the total pre- 
cipitation there. The percentage contribution is still 


highest locally, where on average 3% of the total pre- 
cipitation is from irrigation water during August. 

4. Summary and discussion 

It is well known that irrigation can greatly alter the 
land surface energy and water fluxes, but the contribu- 
tion of irrigation water to precipitation is less well un- 
derstood because of limits in modeling and observations. 
By comparing the evapotranspiration (ET) estimates 
with and without irrigation from the reanalyses them- 
selves and from a hydrological model with detailed 
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a) Annual precip. from irrigation (mm) 




b) Percentage of total precip. (%) 



d) Meridional average (mm/mon) 




0.5 0.75 1 1.25 1.5 1.75 2 2.25 2.5 


Fig. 10. As in Fig. 8, but for irrigation over California. 


considerations of agriculture, this study first demon- 
strates the deficiency of ET estimates from MERRA- 
Land by showing that it has little consideration of the 
impact of irrigation, which can be very large over 
some regions. This enables us to add additional ET 
from irrigation to the ET of MERRA-Land and sepa- 
rate the impact of irrigation on precipitation. The 
quasi-isentropic back-trajectory (Q1BT) method is used 
to determine the evaporative moisture sources for 
precipitation over land using MERRA, but the ET of 
MERRA is corrected using ET from MERRA-Land 
plus estimated ET increase due to irrigation, and 


MERRA precipitation is corrected using CPC unified 
precipitation. 

Results show that the irrigation-caused precipitation 
increase is locally large over some heavily irrigated 
regions, like northern India, eastern Pakistan, and cen- 
tral and northern China, but it is very small in amount 
and percentage for the global average. Irrigation-caused 
precipitation increase is less than the ET increase over 
most areas, indicating that irrigation may be unsustain- 
able over these regions if there are not enough sus- 
tainable water resources. For the same ET increase, 
precipitation increases are greater over wet areas than 
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a) Annual precip. from irrigation (mm) b) Percentage of total precip. (%) 



0.03 0.06 0.09 0.12 0.15 0.18 0.21 0.24 0 . 2 / 0.3 0.33 



0.5 1 1.5 2 

Fig. 11. As in Fig. 8, but for irrigation over the Nile Valley. 


dry areas, but the percentage increases in precipitation 
are similar and are proportional to the ET increase. 
Water vapor from irrigation-caused ET increases over 
arid regions is mostly transported to more humid regions 
where a small increase in humidity can trigger increased 
precipitation. Four heavily irrigated regions are selected 
for detailed analysis. Although they show differences in 
contribution to local and remote precipitation, which 
are related to wind and regional climate, all of them 
show the highest percentage contribution to precipita- 
tion over local land areas. 

As discussed, the ET from MERRA-Land may al- 
ready contain some small contribution from irrigation 


because of the atmospheric data assimilation. This may 
lead to overestimation of the ET without irrigation and 
underestimation of the fraction of ET from irrigation. 
However, the total ET with irrigation is also over- 
estimated. Therefore, the total impact of the bias on the 
estimated irrigation-caused precipitation is uncertain. 
Although there are uncertainties in the method and 
irrigation data and biases in MERRA, this study for 
the first time systematically estimates the direct contri- 
bution of irrigation water to precipitation from a water 
cycle perspective. The estimates in this study are 
calculated under the same atmospheric circulation and, 
so, are complementary to the previous estimates from 
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model simulations and observations, which are not 
suited to separate direct effect from secondary impacts 
and feedbacks. 

This study shows that simple estimates of irrigation 
based on energy limits or soil moisture, which have been 
used by some previous studies, often underestimate or 
introduce seasonal biases in the irrigation impacts on 
ET. This has also been shown by Guimberteau et al. 
(2012; cf. their Fig. 9). If such simple estimates are used 
in climate simulations, the irrigation-caused climate ef- 
fects could be underestimated or incorrectly estimated. 
Our study shows that irrigation can more than double 
ET without irrigation over some regions, which could 
lead to significant climate effects in addition to precip- 
itation changes. Some studies show that the regional 
effect of irrigation could be as large as those seen for 
land cover change (e.g.. Sacks et al. 2009; Douglas et al. 
2009). This suggests that irrigation and its related water 
management should be seriously considered in future 
projections of climate change. 
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